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ABSTRACT 


A  procedure  for  calculating  the  spindle  torque  and  for 

•  » 

determining  the  "optimum"  spindle  axis  location  of  a  controllable- 
pitch  propeller  at  design  conditions  is  presented.  From  the 
blade  geometry  and  theoretical,  pressure  distribution  the  hydro- 
dynamic.  and  centrifugal  forces  acting  on  the  blade  are  obtained, 
and  from  the  distribution  of  these  forces  the  spindle  torque  is 
calculated.  The  optimum  spindle  axis  location  is  then  determined 
by  an  iteration  process.  An  example  ir  included  no  illustrate 
the  procedure. 


•  INTRODUCTION 

* 

A  controllable-pitch  propeller  is  defined  a  propeller  ir  which  the 
pitch  can  be  changed  during  p repel  ler  eperanir-r.  by  turning  earn  blade  about 
a  radial  or  near  radial  axis.  This  pinch- changing  axis  is  called  the 
spindle  axis.  There  are  'wc  ayp-es  of  forces  arcing  on  a  propeller  blade 
under ‘operating  condition;  !I)  hydrodynamic  forces  composed  of  lifting 
forces  (due  to  the  pressure  distribution')  and  viscous  drag  forces,  and 

(2)  centrifugal  forces  (due  to  the  rotating  blade  mass).  The  distribution  of 

•  •  *  • 

these  forces  over  thf • blade  produces  a  ’crque  about  the  spindle  axis.  This 

torque  is  called  ‘he  ,-pi'die  rque  If  10  desirable  tc  have  the  maximum 

* 

spindle  forqj'-  encountered  uncer  •  rainary  ~pera‘ ir.g  conditions  as  small  a; 

•  • 

possible  ir  order  tc  minimise-  "he  maximum  V.sa.-  placed  upon  the  pitch¬ 
changing  mechanism.  This  report,  pre -  an  appre  ximate  method  of  locating 
the  spindle  axis  of  any  given  cladr  lO  as  *.o  produce  zero  spindle  torque 
at  design  condition'-'-  end,  n.-r  cv v r-  uTat  . ••  -r  1*. ch  propellers  with  previously 
located  spindle  axe-,  of  calculating  vh<-  •  pir.dle  t- rque  at  design  conditions 
A  controllable-pitch  propeller  :  '  uatigr.sd  by  essentially  the  same 
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of  pressure  distribution  and  centrifugal  force  distribution  on  a  propeller 
blade  over  the  range  of  operating  conditions  (forward,  backing,  and  tran¬ 
sient  conditions)  follow  different  laws,  the  hydrodynamic  and  centrifugal 
components  of  the  blade  spindle  torque  cannot  be  made  to  balance  over  the 
entire  operating  range  .  However,  the  parameters  can  be  chosen  so  as  to 
minimize  the  largest  spindle  torque  that  will  be  encountered  under  any 
operating  condition  provided  the  effects  of  the  various  parameters  upon 
the  magnitude  of  the  spindle  torque  are  known  throughout  the  range  of 
operating  conditions  and  a  method  of  calculating  the  spindle  torque  for 
any  condition  is  available . 

Additional  compromises,  depending  upon  the  application  of  the  pro¬ 
peller,  must  be  made  in  the  propeller  design.  For  harbor  tugs,  in  which 
backing  is  an  important  consideration,  some  sacrifice  of  forward  performance 
must  be  made  in  order  to  obtain  better  backing  characteristics.  Symmetrical 
or  slightly  cambered  sections  may  be  employed  to  obtain  a  more  favorable 
pressure  distribution  and  better  cavitation  characteristics  under  backing 

conditions.  A  pitch  distribution  with  reduced  pitch  at  inner  sections  of 
*  «* 

the  propeller  blade  at  forward  pitch  settings  is  necessary  to  prevent  the 

inner  blade  sections  from  having  a  positive  pitch  at  negative  nominal 

•pitch  settings.  For  propellers  producing  greater  speed  and  power,  large 

compromises  in  camber  and  pitch  distribution  cannot  be  made  without  great 

sacrifice  in  full  speed  a*head  performance  due  to  cavitation  considerations 

at  full  speed  ahead.  The  inner  sections  of  the  propeller  blade  must  be 

% 

made  thicker  than  for  comparable  fixed-pitch  propellers  due  to  the  unusually 
high  blade  loadings  encountered  under  transient  conditions  during  emergency 
stopping  and  backing  operations. 

In  order  to  permit  the  housing  of  the  pi^ch- changing  mechanism,  the 

hub  size  must  be  somewhat  large  on  a  controllable-pitch  propeller  than 

on  a  fixed-pitch  propeller  (hub  diameter  ratio  of  order  of  magnitude  of 

.30  to  .35  for  a  controllable  pitch  propeller).1  This  larger  hub  size 

will  cause  only  a  slight  less  in  maximum  propeller  efficiency  (or  order 

of  magnitude  of  2  percent  or  less)1  provided*  t  he  hub  and  blades  are 

properly  faired.  Due  to  space  limitations  near  the  hub,  the  number  of 

2 

blades  must  often  tie  1  imited  to  three . 
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References  are  listed  on  page  29. 


There  is  little  quantitative  information  available  on  the  effects 
of  the  various  parameters  upon  spindle  torque.  Rusetskiy  and  Favorskaia"5 
investigated  the  effect  of  blade  outline  and  nominal  pitch  ratio  upon 
the  hydrodynamic  component  of  spindle  torque  at  dead  pull  condition.  They 
developed  a  relationship  whereby  section  angle  of  attack  is  obtained  for 
any  pitch  setting  at  dead  pull  conditions.  The  section  moment  coef¬ 
ficient  is  then  obtain  from  standard  airfoil  data...  and  by  use  of  a 
"coefficient  of  compensation"  (dependent  upon  the  location  of  the  spindle 
axis),  a  representative  spindle  torquie  is  obtained  from  this  moment  coef¬ 
ficient.  Calculations  are  performed  for  three  blade  outlines  and  several 
nogiinal  pitch  ratios  over  a  complete  range  of  blade  pitch  settings. 

When  the  present  investigation  was  almost  completed,  the  author  dis- 

k 

covered  a  method  of  calculating  spindle  torque  devised  by  Van  de  Voorde . 

Van  de  Voorde  presents  a  somewhat  different  approach  to  the  problem  of 
calculating  the  spindle  torque  than  That  of  the  present  paper.  The  dif¬ 
ference  lies  in  the  methods  by  which  the  calculations  are  performed.  Van 
de  Voorde  calculates  the  hydrodynamic  component  of  blade  spindle  torque  by 
use  of  moment  coefficients  of  standard  NACA  airfoil  sections  and  centrifugal 
component  of  blade  spindle  torque  by  direct  integration  throughout  the.  .■ 
entire  blade  mass  of  the  differential  torque  due  to  the  centrifugal  force 

on*  each  element  of  blade  mass »  ^ 

•  « 

In  the  present  report,  *  re-  optimum  spindle  axis  (i  e.,  the  radial 

axis  about  which  the  net.  torque  due  to  the  forces  on  the  blade  is  zero) 
is  determined  by  an  iteration  process:-  in  ,r.ich  me  spindle  torque  is  cal¬ 
culated  about  successive  approximation.-  *■  th°  >'j*i:num  spindle  axis. 

The  hydrodynamic  component,  of  blade  -tind.e  torque  is  calculated  from 

the  theoretical  pressure  distribution  over  fhe  blade  surface.  The  theo- 

5  6 

retical  distribution  over  each  blade  section  is  obtained,  ’  and  •  r.a  dis¬ 
tance  from  any  point  on  the  blade  to  the  initial  spindle  axis  is  calculated 
from  the  blade  geometry.  The  torque  due  to  the  pressure  at  each  of  a  seri-.-s 
of  points  along  each  section  is  then  calculated,  and  this  differential  rorq-.; 


distribution  is  numerically  integrated,  first  along  each  blade  section  and 
then  radially  from  root  to  tip,  yielding  the  total  hydrodynamic  component 
of  spindle  torque. 

A  method  of  calculating  the  component  of  spindle  torque  d he  to  the 
centrifugal  force  distribution  over  each  blade  section  in  terms  of  the 
moments  of  inertia  of  the  blade  sections  is  developed  and  used.  Numerical 
integration  of  these  differential  spindle  torques  is  performed  radially 
from  root  to  tip,  yielding  the  centrifugal  component  of  blade  spindle 
torque . 

' GENERAL  CONSIDERATIONS 

The  shape  of  the  blade  sections  and  the  operating  conditions  of  the 

blade  fix  the  static  pressure  distribution  and  the  centrifugal  force  dis- 

* 

tributions  over  each  section.  If  a  relationship  approximating  these  dis¬ 
tributions  can  be  obtained  from  the  operating  conditions,  then  the  method 
of  approach  to  find  the  spindle  torque  is  to  integrate  the  force  distribu¬ 
tion  over  each  section  to  find  the  resultant  sectional  torque,  and  then 
integrate  these  sectional  spindle  torques  radially  t  find  the  spindle 
torque  of  the  entire  blade. 

The  hydrodynamic  component  of  blade  spipdle  tc ;\r  e  .s  due  to  the 

chordwise  static  pressure  distribution  over  the  blah,  sections.  The 

pressure,  distribution  over  a  given  airfoil  section  is  a  function  of  the 

angle  of  attack  and  the  velocity  of  the  fluid  flowing  over  it.  Thus  for 

a  propeller  blade  section, the  pressure  distribution  is  a  function  of  the 

rotational  velocity  and  the  velocity  cf  advance. 

For  an  airfoil  section, the  ".basic"  pressure  distribution^  that  is, 

6 

the  pressure  distribution  at  ideal  angle  cf  attack,  is  a  function  of  the 

section  shape.  The  location  of  the  center  of  this  pressure  distribution 

is  slightly  forward  of  the  half-chord  point  for  ordinary  aferfoil  sections 

and  is. a  function  of  the  type  of  mean  line  used.  For  the  uniform  load  line 

(a  .»  1.0  mean  line),  this  center  of  pressure  is  very  near  the  l/2  chord 

point  and  for  mean  lines  with  a  calculated  decrease  in  loading  ever  the 

*sar  portion  of  the  section  (such  as  a  -  0.8  mean  line), this  center  of 

v 

pressure  is  closer  to  the  leading  edge  .. ' 

h 


C 
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The  center  of  the  "additional"  pressure  distribution)  that  is, 
pressure  distribution  due  to  angle  of  attack  as  measured  from  ideal  angle 
of  attack,  is  located  at  approximately  25  percent  of  the  chord  from  the 
leading  edge  for  all  sections  The  magnitude  of  this  additional  pres¬ 

sure  distribution  increases  with  increasing  angles  of  attack  and  thus  the 
center  of  pressure  of  the  sections  moves  towards  the  quarter-chord  point. 

If  the  propeller  is  designed  so  that  the  sections  operate  at  ideal 
angle  of  attack  at  design  operation  conditions,  the  pressure  distribution 
over  each  section  at  design  conditions  will  be  the  basic  pressure  distri¬ 
bution.  As  the  advance  coefficient  decreases,  the  angle  of  attack  becomes 
larger  than  design  angle  of  attack,  the  additional  pressure  distribution 
increases  in  magnitude,  and  thus  the  center  of  pressure  on  the  blade 
sections  moves  toward  the  quarter-chord  point. 

In  the  design  of  a  propeller,  the  desired  lift  and  hydrodynamic 

pitch  angle  at  each  radius  are  determined  and  then  a  section  is -selected 

to  produce  the  desired  lift  in  two-dimensional  flow.  However,  if  the 

propeller  were  built  with  the  two-dimensional  sections  (same  camber  and 

angle  of  attack)  laid  along  the  helices,  the  propeller  would  be  found  from 

experience  to  be  underpitched-  The  deficiency  in  pitch  is  due  to  "three- 

dimensional  effects"  in  the  propeller.  These  three-dimensional 'effects 

8 

are  compensated  for  in  the  design  by  applying  correction  factors  to  the 
camber  and  angle  of  afack  a-  each  section.  With  these  correction 
factors,  propellers  are  designed  that  develop  the  desired  overall 
characteristics  (thru.*,  pitch,  cavita'ion  suppression,  etc.),  but  the 
lift  distribution  over  each  section  is  net  rec“i?arily  the  same  as  that 
over  the  two-dimensional  section  from  which  it  was  developed. 

The  actual  lift  distribution  approximates  the  two -dimeric ions  1  life? 
distribution  only  to  the  degree  that  the  correction  factors  adjust  the 
two-dimensional  sections  to  the  actual  three-dimensional  flow  r*' stern 
along  the  blade  section.  In  -'re  present  propel, .er  design  met  nod,  the 
corrections  are  based  only  •  n  certain  positions  * 1  the  section,  no* 
continuously  along  the  section-  Hence,  -here  is  no  way  ci  theoretically 
determining  the  true  three- iimerv  tonal  lift  distribution,  and  since  no 
accurate  experimental  measurement  of  the  actual  pressure  distribution 


* 


existing  over  a  propeller  blade  has  yet  been  made,  it  will  be  assumed  that 
the  lift  distribution  is  to  be  the  same  as  the  two-dimensional  lift  dis¬ 
tribution;  that  is,  the  correction  factors  applied  in  the  propeller  design 
are  assumed  to  exactly  adjust  the  section  from  the  two-dimensional  flow  to 
three-dimensional  flow . 

In  the  calculation  of  the  pressure  distribution  the  blade  sections 
are  assumed  to  operate  continually  at  the  mean  inflow  velocity  and  angle 
of  attack.  But  in  actuality  there  is  a  wake  variation  which  might,  in 
the  course  of  one  revolution  of  the  propeller,  cause  significant  variation 
in  the  sectional  angles  of  attack  and  hence  in  the  pressure  distribution 
and  spindle  torque .  This  could  cause  instantaneous  spindle  torques  of  much 
greater  magnitude  than  those  calculated  at  the  mean  flow  conditions .  Thus 
a'  reasonable  factor  of  safety  should  be  applied  to  any  spindle  torque 
calculations  based  upon  mean  flow  conditions. 

The  centrifugal  component  of  spindle  torque  arises  from  the  centrifugal 
force  acting  radially  outward  on  each  increment  of  the  blade  mass  of  the 
rotating  propeller.  This  component  of  spindle  torque  is  due  to  the  fact 
that, in  general,  the  centrifugal  force  on  each  increment  of  blade  mass 
has  a  component  in  the  plane  normal  to  the  spindle  axis. 

Such  design  factors  as  the  amount  of  blade  rake  and  skew  have  a  large 
effect  upon  the  magnitude  of  the  blade  spindle  torque  as  they  offset  the 
sections  from  the  radial  disc  line .  Large  amounts  of  rake  and  skew  are 
not  generally  used  for  controllable-pitch  propellers  due  to  clearance 
limitations  over  the  range  of  pitch  settings. 

The  spindle  axis  is  usually  radial,  but  for  applications  in  which 
stern  clearance  is  a  problem  at  negative  pitch  settings  the  spindle  axis 
is  Sometimes  sloped  in  the  af^ward  direction.  For  the  purposes  of  cal¬ 
culations  in  this  report  the  spindle  axis  i ■  considered  +o  be  radial. 

*  "  PROCEDURE 

The  primary  problem  considered  is;  given  a  made  desigft  without  the 
spindle  'located,  find  the  location  of  the  spindle  so  that  the  torque  about 
the  spindle  axis  due  tc.  the  forxs  acting  the  blade  at  design  conditions 
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is  zero.  This  "optimum"  spindle  axis  location  is  determined  by  an  itera¬ 
tion  process.  The  torque  is  first  calculated  at  mt  a  selected  abstract 
axis  and  from  the  magnitude  of  this  torque,  a  second  abstract  axis  is 
determined  which  better  approximates  the  optimum  spindle  axis  location. 

S1.  •'cessively  closer  approximations  to  the.  optimum  spindle  axis  location 
are  made  until  the  calculated  torque  about  an  axis  is  sufficiently  near 
zero,  this  final  axis  being  the  determined  optimum  spindle  axis  location. 
Thus  to  obtain  zero  torque  at  design  conditions  the  spindle  should  be 
located  so  that  the  actual  axis  of  twist  is  this  optimum  spindle  axis. 

As  a  convenient  starting  point,  the  first  abstract  axis,  called  the 
"initial  spindle  axi^"  is  taken  as  the  projection  of  the  pitch  reference 
line  in  a  plane  normal  to  the  propeller  shaft  and  intersecting  the  pitch 
reference  line  in  the  section  at  the  hub.  This  location  is  selected 
because  the  pitch  reference  line  passes  through  the  chord  line  of  each 

blade  section  And  serves  as  the  reference  line  when  delineating  a  screw 

9 

propeller.  Thus  a  relatively  simple  expression  for  the  distance  to  this 

initial  spindle  axis  from  any  point  on  the  blade  can  readily  be  obtained 

from  blade  geometry.  For  subsequent  abstract  spindle  axis  locations, 

the  expressions  for  distance  are  slightly  more  involved.  For  unskewed 

blades  with  maximum  thickness  at  the  section  half-chord  point,  the  pitch 

reference  line  passes  through  the  half-chord  point  of  each  section;  for 

skewed  blades  it  usually,  but  not  always,  passes  through  the  half-chord 

9 

point  of  the  section  at  the  hub. 

A  second  type  problem  that  can  also  be  solved  by  the  method  of  this 
report  is:  given  a  completely  designed  blade  with  the  actual  spindle  axis 
location  already  fixed,  find  the  spindle  torque  at  design  conditions.  In 
this  case,  the  spindle  torque  is  directly  calculated  about  the  actual 
spindle  axis . 

HYDRODYNAMIC  COMPONENT  OF  BLADE  SPINDLE  TORQUE 

The  hydrodynamic  component  of  spindle  torque  is  calculated  from  the 
pressure  distribution  over  the  blade  surface .  The  pressure  distribution 
may  be  determined  by  any  available  means  since  it  is  an  input  to  the 
method  of  spindle  torque  calculation  .  For  NACA  sections  the  pressure 
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Figure  I  -  Propeller  Isometric  Showing  Hydrodynamic  Forces 
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SECTION  AT  HUB 


Figure  2  -  Diagram  of  Hydrodynamic  Forces 
and  Distances  to  Initial  Spindle  Axis 


distribution  may  be  obtained-  quickly  by  the  MCA  approximate  method. 

For  such  sections  operating  at  ideal  angle  of  attack,  the  pressure  dis¬ 
tribution  and  center  of  pressure  are  functions  only  of  the  type  mean  line 
and  may  be  obtained  directly  from  MCA  data.  For  propellers  with  blade 

sections  which  are  not  of  MCA  type,  the  pressure  distribution  may  be 

5 

obtained  by  the  AVA  approximate  method.  This  method  may  be  used  to 

calculate  the  pressure  distribution  over  a  body  of  arbitrary  shape  from 

the  ordinates  measured  from  the  longest  chord  contained  within  the  body 

and  the  operating  angle  of  attack  measured  from  this  chord.  The  big 

disadvantage  to  this  method  is  the  tediousness  of  the  work  involved. 

Although  the  calculations  are  programmed  for  the  TJNIYAC  high-speed 
10 

computer,  the  ordinates  must  be  measured  on  large-scale  drawings 
(48"  chord)  for  nonmathematieai  sections,  and  the  pressure  is  obtained  at 
inconvenient  stations  along  the  chord. 

The  lift  and  center  of  pressure  of  each  blade  section  are  obtained 
by  integration  of  the  pressure  distribution  (assumed  to  act  normal  to  the 
velocity  past  the  section,  this  velocity  being  the  resultant  inflow 
velocity  to  the  section).  The  viscous  drag  coefficient  is  obtained 
from  MCA  dataj^  and  the  drag  force  acts  at  the  center  of  pressure  in  the 
direction  of  the  velocity  past  the  section.  These  lift  and  drag  forces 
for  each  blade  section  are  tbr-r.  decomposed  into  axial  and  tangential  com¬ 
ponents.  From  the  geometric  characteristics  of  the  blade,  expressions  are 
derived  for  the  distances  from  the  initial  op  inale  axis  to  any  axial  and 


•’T’no  e} *■*•'*-  r"-  ■?  2  W 


due 


tangential  force  on  the  blad<- 
lift  and  drag  on  each  blade  section  is  then  computed.  These  differential 
spindle  torques  are  then  -radially  integrated  ever  the  blade  from  root  to 
tip,  yielding  the  hydrodynamic  component  of  spindle  torque  about  the 
spindle  axis.  The  differential  hydrodynamic  component  of  blade  spindle 
torque  at  each  radius  is  given  by  (see  Figures  1  and  2): 


Qf,  fa*  (^l  -  Cj,  SM/i{)i-V(cL  s,n/6<'  tCpcos/Scj]^  V*/  dr 

* 


[lj 


where  u  and  v,  the  tangential  and  axial  components,  respectively,  of  the 
distance  from  the  initial  spindle  axis  to  the  center  of  pressure  are  given 
by: 


For  any  given  two-dimensional  pressure  distribution,  the  lift  coefficient 


Csfirf dc. 


and  the  center  of  pressure  (in  percent  chord  from  leading  edge)  is 


cc-r. 


cxjcx 


The  drag  coefficient,  C^,  for  standard  airfoil  sections  may  be  obtained 
from  experimental  NACA  data.  At.  design  conditions,  the  drag  contribution 
to  spindle  torque  will  be  very  small,  since  the  blade  sections  are  generally 
designed  to  operate  at  ideal  angle  of  attack  where  the  drag  coefficient 
is  a  minimum.  . 

At  off  design  conditions,  the  sections  will  net  in  general  operate  at 
design  angle  of  attack,  subsequently  the  drag  coefficients  may  be  many 
•times  greater  tha.i  at  design  conditions  and  the  spindle  torque  due  to  drag 
may  become  more  significant . 

The  dynamic  pressure,  q,  -is  based  upon  the  resultant  inflow  velocity, 

V  ,  hence: 

T* 
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where 


v;  =14 


CQSfft;  • 

sm/S 


for  free  running  propellers 


Vr*VO-v,) 


C0S(/3t- 

sisifi 


for  wake  adapted  propellers 


From  the  propeller  design  calculations,  the  values  of  V  w  ,  /S^,  ft, 
and  o<  are  known  at  design  conditions.  The  only  other  needed  quantities 
are  geometric  characteristics  of  the  propeller  (R  ,  S^,  1,  <f>  ,  C  ,  R) 
which  can  easily  be  obtained  from  the  propeller  drawing  or  design  calcu¬ 
lations.  The  sign  conventions  are  chosen  so  that  spindle  torques  tending 
to  rotate  the  blade  toward  larger  forward  pitch  settings  are  positive . 

R  is  the  offset  of  a  section  from  the  spindle  axis  due  to  rake  and  is  thus 

9. 

the  offset  of  the  section  from  the  propeller  disc  minus  the  offset  of  the 
section  at  the  hub  from  the  propeller  disc.  Hence  for  a  propeller  with  a 
rake  angle  ^  ,  the  offset  of  the  section  at  radius  due  to  rake  is  simply 

=  (  r- rh)ta.n  J  [8] 

The  hydrodynamic  component  of  blade  spindle  torque  is  then  obtained 
by  integrating  Equation  1  radially  over  the  blade  from  root  to  tip: 

Qh  ~  /~RJ  Lv(cl  cos  ft ^  -Gj  si  *y3t)  +  y(CL  -tQ  cos^jj  Vr  Jx  [  9  ] 


CENTRIFUGAL  COMPONENT  OF  BLADE  SPINDLE  TORQUE 

An  expression  is  derived  which  expresses  the  differential  centrifugal 
component  of  spindle  torque  at  each  radius  In  terms  of  the  maximum  and 
minimum  moments  of  inertia  of  the  blade  sections,  the  blade  sectional  area, 
and  the?  axial  and  tangential  component  distances  from  the  spindle  axis  to 
the  blade  section  centroid.  The  maximum  and  minimum  sectional  moments  of 
inertia,  the  location  of  section  centroids,  and  the  sectional  areas  are 
estimated  by  the  mtethcd  in  Reference  11.  From  the  geometric  characteristics 
of  the  blade,  expressions  are  derived  for  the  axial  and  tangential  component 
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INITIAL  SPINDLE  AXIS 


Figure  3  -  Propeller  Isometric  Showing  Centrifugal  Forces 


distances  from  the  spindle  axis  to  the  section  centroids.  Radial  integra¬ 
tion  of  these  differential  spindle  torques  yields  the  total  centrifugal 
component  of  blade  spindle  torque. 

The  differential  centrifugal  component  of  spindle  torque  at  each 
radius  is  approximated  by; 


+  FIY 1 


J 


a  * 

r 


[10]  - 


The  values  of  Y  and  Z  in  Equation  10  relative  to  the  initial  spindle  axis 
are  calculated  from  the  geometric  characteristics  of  the  blade  as  follows; 


Y-  r  s/m 


To4^)a  ~  sk  YXcj  coscf)  -Yc  SIK ij) 


[n] 


+  Yc  cos  <j)  - 


[12] 


Reference  11  presents  a  simplified  approximate  method  for  calculating 
I  ,I.,X,Y,  and  A  for  airfoil  sections  of  the  type  commonly  used 
in  propeller  design.  For  other  type  sections  these  values  may  be  obtained 
by  numerical  integration.  These  sectional  geometric  characteristics  in 
Equations  10,  11,  and  12  refer  to  the  projection  of  the  developed  sections 
in  a  plane  normal  to  the  spindle  axis,  but  the  values  for  the  developed 
sections  may  be  used  as  close  approximations.  The  other  geometric  charac¬ 
teristics  in  Equations  10,  11,  and  12  can  readily  be  obtained  from  the  pro¬ 
peller  drawing.  The  total  centrifugal  blade  spindle  torque  is  then  obtained 
by  integration  of  Equation  10  over  the  blade  from  hub  to  tip; 

[13] 

/'h 

*  The  derivation  is  presented  in  the  appendix. 
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Figure  5  -  Effect  of  Changed  Sp 
on  Hydrodynamic  Blade  Sp 
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PLAN  VIEW 


Figure  6  -  Effect  of  Changed  Spindle  Axis  Location 
on  Centrifugal  Blade  Spindle  Torque 
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CHANGED  LOCATION  OF  SPINDLE  AXIS 


After  the  spindle  torque  has  been  calculated  about  the  initial 
spindle  axis,  successive  spindle  axis  locations  are  chosen  and  the  spindle 
torque  calculated  about  these  axes  in  an  effort  to  find  the  optimum  spindle 
axis  location.  For  convenience  of  calculation,  these  axes  are  taken  radi¬ 
ally  through  the  chord  of  the  section  at  the  hub  in  a  plane-  normal  to--the 
propeller  axis.  Thus  all  axes  lie  in  the  helical  surface  with  pitch 
equal  to  that  of  the  section  at  the  hub  at  design  pitch  setting.  Hence 
each  axis  can  be  regarded  as  rotated  from  the  initial  spindle  axis  by  an 
angle,  say  © ,  about  the  propeller  axis  and  shifted  axially  from  the 
initial  spindle  axis,  a  distance  cf  r^0 tan^i^. 

The  distances  from  the  new  spindle  axis  location  to  the  components 
of  forces  on  the  blade  are  then? 

* 

(l)  for  the  hydrodynamic  blade  spindle  torque 


v«  SIN<P-^-ri i©  +atlh  [151 


(2)  for  the  centrifugal  blade  spindle  torque 


r 
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2  SK  tXc]siN(j:  +  vc  Cr5(f  ~Z.  -  rh  a  to.  [l7j 

Shifting  the  spindle  axis  towards  the  leading  edge  decreases  the 
distance  to  the  lift  forces  near  the  leading  edge  of  the  sections  and 
increases  the  distance  to  the  j-iff  forces  near  the  trailing  edge  of  the 
sections.  This  has  the  effect  of  shifting  the  hydrodynamic  spindle  torque 
at  design  conditions  toward  negative  values  (tending  towards  lower'  posi¬ 
tive  pitch  settings/. 


Since  the  direction  of  movement  of  the  spindle  axis  is  in  a  helical 
surface  with  pitch  equal  to  that  of  the  hub  section,  the  spindle  axis 
movement  at  each  radius  is  in  the  general  direction  parallel  to  the  section 
chord  at  that  radius  (see  Figure  5 ) °  hue  to  the  radial  pitch  variation, 
it  is  not  actually  parallel  to  the  chord.  The  direction  of  spindle  axis 
movement  is  nearly  perpendicular  to  the  lift  forces,  and  thus  movement 
is  in  the  direction  so  as  to  have  the  maximum  effect  on  the  hydrodynamic 
blade  spindle  torque. 

The  centrifugal  blade  spindle  torque  is'  due  to  the  component  of 
centrifugal  force  in  the  plane  normal  to  the  spindle  axis.  When  the  spindle 
axis  is  moved,  the  distance  to  each  increment  of  centrifugal  force  is 
changed  by  the  axial  distance  which  the  spindle  axis  is  moved  (see  Figure 
6)  and  the  components  of  centrifugal  force  in  the  plane  normal  to  the 
spindle  axis  are  slightly  changed  due  to  the  angular  movement  of  the 
spindle  axis .  The  centrifugal  component  of  spindle  torque  is  thus  expected 
to  be  much  less  sensitive  to  spindle  axis  location  than  is  the  hydrodynamic 
component  of  blade  spindle  torque. 

Shifting  the  spindle  axis  towards  the  leading  edge  of  the  section 
at  the  hub  at  design  pitch  setting  shifts  both  components  of  the  section 
centroid  offset  from  the  spindle  axis  (Y  and'Z)  towards  negative  values. 

For  blades  with  after  rake  and  positive  skevback,  both  Y  and  Z  are  nega¬ 
tive  and  shifting  the  spindle  axis  towards  the  leading  edge  increases  the 
negative  centrifugal  component  of  spinal*-  torque  due  fc  the  offset 
centroid . 

Thus,  for  a  propeller  with  after  rake  ard  positive  skevback, shifting 
the  spindle  axis  towards  the  leasing  edge  of  the  hub  section  will  cause 
both  the  hydrodynamic  and  centrifugal  component-  of  spindle  torque  to 
change  toward  negative  values  and  the  change  in  the  hydrodynamic  component 
can  be  expected  to  be  much  greater  than  the  change  in  centrifugal  component. 

OFF-DESIGN  CONDITIONS 

* 

The  method  outlined  ir.  this  report,  for  calculating  the  spindle  torque 
*  *  *  . 

may  be  applied  over  the  complete  rang*.-  of  operating  conditions  (range  of 


loadings  and  pitch  settings)  provided  that  the  pressure  distribution  and 
centrifugal  force  distribution  are  known  at  each  condition.  By  the  cal¬ 
culation  of  spindle  torque  about  several  spindle  axis  locations  over  the 
complete  range  of  operating  conditions,  the  spindle  may  be  located  so  that 
the  maximum  spindle  torque  encountered  under  steady-state  conditions  will 
be  as  small  as  possible,  regardless  of  the  operating  condition  at  which 
this  maximum  occurs . 

The  pressure  distribution  over  a  blade  section  is  a  function  of  the 
effective  geometry  of  the  section  (thickness  and  camber)  and  the  angle  of 
attack  at  which  the  section  is  operating.  At  a  given  pitch  setting,  the 
geometry  of  each  blade  section  remains  unchanged  but  the  angle  of  attack 
at  which  each  section  is  operating  varies  as  a  function  of  the  propeller 
advance.  References  12  and  13  give  approximate  methods  of  calculating  the 
section  angles  of  attack  at  off- design  advance.  These  methods  are  meant 
to  apply  to  fixed  pitch  propellers  but,  for  lack  of  a  better  method,  they 
could  be  used  for  a  controllable-pitch  propeller  at  off-design  pitch 
settings.  As  the  blade  is  rotated  about  the  spindle  axis,  the  section  of 
the  blade  cut  by  a  cylindrical  surface,  whose  axis  coincides  with  the 
propeller  axis,  changes.  Since  for  flow  considerations,  a  blade  section  is 
considered  as  located  at  a  constant  radius,  the  effective  blade  sections 
change  as  the  blade  is  rotated  about  the  spindle  axis.  This  complicates 
the  problem  of  determining  pressure  distribution  at  off-design  pitch 
settings.  For  a  helical  surface  rotated  about  a  radial  line.  Reference 
lU  gives  a  method  of  determining  the  deformation  of  this  surface  from  a 
helical  surface  through  two  reference  points  of  the  original  surface  in 
its  displaced  position,  at  each  radius  in  question. 

The  centrifugal  component  of  spindle  torque  is  a  function  of  pitch 
setting  and  angular  velocity  but  is  not  a  function  of  blade  loading.  The 
effective  change  in  blade  section  shape  is  small  and  its  effect  upon  the 
centrifugal  spindle  torque  is  probably  insignificant.  If  the  effect  of 
variation  in  section  geometry  is  ignored,  the  centrifugal  component  of 
spindle  torque  varies  only  as  a  function  of  the  propeller  angular  velocity 
and  the  pitch  angle  of  each  blade  section:  thus  it  can  be  calculated  over 
the  complete  range  of  operating  conditions  without  great  difficulty. 


\ 
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For  applications  in  which  the  pitch  setting  must  be  reversed 
rapidly  for  maneuvering  purposes,  the  magnitude  of  the  spindle  torque 
during  transient  periods  becomes  important.  Due  to  the  complex  and 
unsteady  nature  of  the  flow  during  emergency  ahead  and  emergency  astern 
operations,  it  is  difficult  to  estimate  the  pressure  distributions  for 
these  conditions .  Experimental  study  seems  to  be  the  best  approach  to 
this  problem. 


EXAMPIE 

An  example  calculation  will  now  be  performed  to  illustrate  the 
procedure . 

The  calculations  will  be  performed  on  the  propeller  from  Appendix 
2  of  Reference  8  in  which  the  design  calculations  for  the  propeller 
are  performed.  By  referring  to  Reference  8,  the  reader  can  more  directly 
see  what  quantities  are  available  for  the  calculation  of  the  spindle 
torque  and  how  these  quantities  are  obtained.  This  example  propeller  is  a 
fixed-pitch  propeller  but  at  design  conditions  the  geometry  of  a  controllable- 
pitch  propeller  is  almost  identical  to  that  of  a  fixed-pitch  propeller 
designed  to  operate  at  the  same  conditions. 

The  design  conditions  for  this  propeller  are  as  follows: 

V  =21  knots 

Pg  =  13,000  ehp 
=  17,500  shp 

rpm  =  102 

to  =10 .68  rad/  sec 
Diameter  =  21  ft 

w  =  0.20 
o 

t  =  0.15 

The  type  blade  section  used  in  this  propeller  is  the  NACA  66  section 
with  parabolic  tail  and  a  =  0.8  mean  line. 

The  problem  considered  in  this  example  is:  given  the  fully  designed 
propeller,  find  the  location  of  the  spindle  so  that  the  spindle  torque 
will  be  a  minimum  at  design  conditions.  As  a  good  starting  point,  tb® 
torque  will  be  calculated  about  the  initial  spindle  axis. 
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A)  HYDRODYNAMIC  BLADE  SPINDLE  TORQUE 

Properties  of  the  additional  blade  sections  at  x  =  .925*  «95> 
and  .975  were  calculated  in  order  to  get  a  clearer  picture  of  the 
radial  spindle  torque  distribution  in  the  region  near  the  blade  tip. 

Step  1  -  The  center  of  pressure  of  the  blade  sections  is  determined. 
This  propeller  design  is  based  upon  shock  free  entry  so  that  all  the  lift 
is  produced  by  camber.  Hence,  assuming  the  propeller  sections  to  operate 
at  the  conditions  'for  which  they  are  designed,  the  pressure  distribution 
over  each  section  is  simply  that  produced  by  the  NACA  0.8  mean  line  as 
given  in  the  NACA  data.  Thus  Cc  =  45.2  for  each  section  (0.8  mean 
line  distribution  and  Equations  4  and  5). 

Step  2  -  The  distances  in  the  axial  and  tangential  directions  from 
the  initial  spindle  axis  to  the  section  center  of  pressure  are  determined 
( columns  6  and  7 ) » 

Column  (See  page  23) 

2  1  from  design  calculation 

3  from  design  calculations 

4  ^  from  design  calculations 

5  R  Equation  (8);  from  design, <5  =  7*5° 

a. 

6  u  Equation  (2);  C  =  45  for  NACA  66  section 

7  v  Equation  (3) 

Step  3  -  The  section  hydrodynamic  component  of  blade  spindle  torque 
is  then  calculated 


(l-v^)  from  design  calculations 
5 from  design  calculations 
siryS  from  design  calculations 

Equat ion  ( 7b } 

C^  from  design  calculations 

cosyd^  from  design  calculations 

sinyS.  from  design  calculations 

Q’h  Equation  1,  -■  .008,  p- 


lb-sec^ 


Step  4  -  The  hydrodynamic  component-  if  blade-  spindle  t-orque  is 
obtained  by  integration  of  Equation  1  from  root  to  tip. 

16  Simpson's  Multiplier 

17  Simpson's  Product  -  column  ip  times  column  l6 

iq  (  2_  coiumn  .  T ) 


b)  centrifugal  blade  spindle  torque 

The  NACA  66  section  with  parabolic  tail  and  a  =  0.8  mean  line  is  used 
in  this  propeller.  The  equations  for  estimating  the  geometric  character¬ 
istics  of  this  section  are  given  on  page  358  of  Reference  8.  The  following 
quantities  are  estimated  by  these  equations: 

A’  =  .963  (at  / l)  1^  (a  from  Figure  24,  Reference  8)  [l8] 

x  -  J/2-  x* 
c,  1 

x  =  (.500  -  A73  +  .026  m  /l)  1  [19] 

Yc  “  "yl 

=  +(.155  t  /l  +  0.80)  (m  /l)  1  .  [20] 

I  .  =  I  (assuming  the  major  axis  of  the  section  is  parallel  to  the 

nose-tail  line] 

=  .9^5  b  (mx/l)2  +  .04487  (t^/l)3  l4  •  [21] 

(b  from  Figure  25,  Reference  8) 

I  =  I  (assuming  the  minor  axis  of  the  section  is  perpendicular 
max  yo  to  the  nose-tall  line) 

=  .914  (c)  (tx/l)  1^  (c,  from  Figure  26,  Reference  8)  [22] 

The  geometric  properties  of  each  blade  section  are  then  calculated 

as  follows: 

Column  (See  page  25) 

2  t x /l  known  from  design  calculations 

3  m  /l  known  from  design  calculations 

X  • 

4  1*  known  from  design  calculations 

5  a  Figure  24,  Reference  $  * 

6  b  Figure  25,  Reference  8 

7  c  Figure  26,  Reference  8 

8  A  Equation  18 

9  Equation  19 

10  Yc  Equation  20 

11  I  .  Equation  21 

mm 

12  I  Equation  22 

max 

The  location  of  the  section  centroid  relative  to  the  spindle  axis  is  then 
found . 

i 13  Y  Equation  II 

14  Z  Equation  12 

24 


I 


1  2  S 

3  4 

5 

6 

7 

8 

9 

10 

11 

x  tx/l  mx/l  1 

a 

b 

c 

A 

X 

c 

Y 

c 

I  ft) 

(ft2) 

(ft) 

(ft) 

(ft  ) 

0.2  0.1900  0.0349  4.620 

0.7481 

0.0097 

0.04355 

2.922 

0.1294 

0.1325 

0.1?.?5 

0.3  0.1434  O.O36O  5-3229 

0-7484 

O.OO65 

0.04350 

2.826 

0.1464 

0.1537 

0.0994 

0.4  0.1111  0.0367’  5.691 

■>  «* 

0.7485 

0.0046 

0.04348 

2.635 

0.1593 

0.1698 

0.0692 

0.5  0.0873  0.0370  5.964 

0.7485 

0.0036 

0.04348 

2.239 

0.1670 

0.1787 

0.0416 

0.6  0.0688  0.0363  6.069 

0.7485 

0.0026 

0.04550 

1.827 

0.1699 

0.1780 

0.0231 

0.7  0.0539  0.0345  5.964 

0.7480 

0.0022 

0.04355 

1.382 

0.1670 

O.1658 

0.0115 

0.8  0.0423  .  0.0320  5-460 

0.7478 

0.0016 

0.04365 

0.908 

0.1526 

0.l4o6 

0.0042 

0.9  0.0330  0.0291  4.389 

0.7475 

0.0013 

0.04375 

0.458 

0.1229 

0.1026 

0.0010 

1.0 

• 

12  13 

14 

15 

16 

• 

17 

18 

19 

*max  1 

Z  ^max 

'  Wsin  2(f) 

AYZ 

V 

SIMPSON’S 

MULTIPLIERS 

SIMPSON'S 

PRODUCTS 

2 

c 

(ft4)  (ft). 

(ft) 

(ft*) 

(ft4) 

(ft-lbs/ft) 

3.444  -0.037 

♦0.188 

♦1.502 

-0.0203 

-2759 

1 

-2759 

4.262  +0.103 

♦0.089 

+2.070 

+0.0026 

-3860 

4 

-15,440 

4.778  *0.035 

-0.022 

+2.328 

-0.0020 

-4330 

2 

-8660 

4.389  *0.008 

-0.190 

♦2.050 

-0.0036 

-3809 

4 

-15,236 

3.709  -0.091 

-0.400 

♦1.640 

♦0.0662 

-3176 

2 

-6552 

2-714  -0.282 

-0.652 

♦1.106 

♦0.2536 

-2529 

4 

-10, lie 

1.581  -0.594 

-0.948 

♦0.598 

♦0.5120 

-2066 

2 

-4132 

0.489  -1.064 

-1.294 

♦0.344 

+0.6304 

-1813 

4 

-7252 

0 

1 

0 

-69,9l" 

%  =  9,947)  =  24,481  ft-ibs. 
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The 

15 

16 
IT 


spindle  torque  at  each  hlade  section  is  now  computed 

(W_~Jmin)  gin  2(f) 

A  YZ 


Q'  -  Equation  10. 

'  A  *  *  Ss 


The  blade  material  is  naval  brass 


The  distribution  of  centrifugal  spindle  torque  is  next  integrated  over  the 
blade  from. hub  to  tip  to  find  the  total  centrifugal  blade  spindle  torque. 

18  Simpson's  Multiplier 

19  Simpson's  Product-  column  17  times  column  l8 


C)  LOCATION  OF  THE  pPTIMUM  SPINDLE  AXIS 

An  approximation  to  the  optimum  spindle  axis  location  is  now  made 
based  upon  blade  geometry  and  the  magnitude  of  the  spindle  torque  about 
the  initial  spindle  axis.  The  spindle  torque  about  this  new  axis  location 
is  then  calculated  and,  along  with  the  previously  calculated  torque  about 
the  initial  spindle  axis  location,  is  used  to  obtain  a  better  estimation 
to  the  optimum  spindle  axis  location.  Successive  calculations  are  per¬ 
formed  until  a  sufficiently  close  approximation  +o  the  optimum  location 
is  obtained;  that  is,  until  the  calculated  spindle  torque  about  an  axis 
is  sufficiently  near  zero. 

In  the  example  calculation,  the  torques  about  the  arbitrary  spindle 
axis  are  found  to  be: 

^  =  -56,146  ft- lb 
qc  =  -24,481  ft -lb 
Qb  +  Qc  =  -60,627  ft-lb 

Thus  in  an  effort  to  find  the  optimum  spindle  axis  location,  the  spindle 
axis  must  be  shifted  somewhat  towards  the  rear  of  the  blade .  As  an- 
approximation  to  the  optimum  spindle  axis  location,  the  second  spindle 
axis  is  taken  radially  through  the  hub  section  at  an  angle  of  4  degrees 
towards  the  trailing  edge,  as  measured  from  the  initial  spindle  axis. 
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Nond i mens i ona  I  Radius,  x 

Figure  9  -  Radial  Distribution  of  Hydrodynamic  and  Centrifugal 
Components  of  Blade  Spindle  Torque  for  Example  Propeller 


Figure  10  -  Variation  of  Spindle  Torque  with  Location 
of  Spindle  Axis  for  Example  Propeller 
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This  is  a  shift  from  the  initial  spindle  axis  of  5*9  percent  of  the 
*• 

huh  section  length. 

The  spindle  torque  about  the  new  spindle  axis  is  calculated  by  the 
same  procedure  as  already  shown  for  the  previous  spindle  axis  location. 

The  only  values  that  have  changed  are  the  distances  to  the  spindle  axis 
from  various  positions  on  the  blade.  For  hydrodynamic  blade  spindle 
torque,  u  and  v  are  calculated  by  Equations  it  and  15  respectively,  and 
for  centrifugal  blade  spindle  torque,  Y  and  Z  are  calculated  by  Equations 
16  and  17,  respectively. 

The  components  of  spindle  torque  about  this  second  spindle  axis 
location  are  calculated  to  be:  * 

^  =  +  1675  ft-lb 
Qc  =  -.22,433  ft-lb 
%  +  Qc  =  -20,758  ft-lb 

This  smaller  spindle  torque  of  the  same  sign  indicates  that  the 
shift  is* in  the  proper  direction  but  is  not  great  enough.  A  third  spindle 
axis  location  is  then  taken  at  -6.4  degrees  (9-4  percent  of  the  hub  section 

*S 

length)  from  the  initial  spindle  axis .  The  components  of  spindle  torque 
about  this  third  spindle  axis  location  are  calculated  to  be: 

=  27,225  ft-lb 
Qc  =  -23,*317  ft-lb 
+  Qc  =  +  3,908  ft-lb 

The  total  torque  about  this  spindle  axis  is  of  reversed  sign,  hence 
this  location  is  a  little  too  far  aft  along  the  blade.  By  plotting  spindle 
torque  versus  spindle  axis  location  (see  Figure  10),  the  optimum  spindle 
axis  location  (zero  spindle  torque*  is  seen  to  be  6.0  degrees  (8.8  percent 
hub  section)aft  the  initial  spindle  axis. 

CONCLUDING  REMARKS 

This  report  has  outlined  a  workable  method  for  calculating  the  spindle 
torque  and  for  (determining  the  optimum  spindle  axis  location  of  a  con¬ 
trollable  pitch  propeller  at  design  conditions  based  upon  two-dimensional 
chordwise  pressure  distributions. 
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Such  design  parameters  as  rake,  skew,  blade  outline,  type  of  blade 
sections,  and  pitch  distribution  are  seen  to  influence 'the  magnitude  of 
the  spindle  torque .  The  results  of  the  example  indicate  that  the  hydro- 
dynamic  component  of  blade  spindle  torque  is  much  more  sensitive  to 
spindle  axis  location  than  is  the  centrifugal  component. 

The  method  of  calculating  the  spindle  torque  outlined  in  this  report 
can  also  be  applied  to  off- design  conditions  provided  that  the  pressure 
distribution  at  these  conditions  is  known.  The  next-  logical  step  is  to 
develop  a  method  of  approximating  the  pressure  distribution  at  off- design 
conditions  and  by  use  of  the  method  presented  in  this  report,  to  investigate 
the  effect  of  rake,  skew,  blade  outline,  type  blade  sections,  and  pitch 
distribution  on  spindle  torque  over  a  complete  range  of  operating  conditions. 
A  subsequent  report  investigating  the  effect  of  these  parameters  over  a 
range  of  operating  conditions  is  planned.  It  is  felt  that  such  an  investi¬ 
gation  would  greatly  facilitate  the  design  of  controllable -pitch  propellers 
*for  minimized  spindle  torque. 


REFERENCES 

1.  Rupp,  L„A„,  "Contra  1 lable-Prtch  Propellers,”  Society  of  NaVai 
Architects  and  Marine  Engineers,  VoL .  56,  pp, 272-358  ( 19^8 ) 

,  2.  Taylor,  J„LV  "The  Variable-Pitch  Marine  Propeller,"  The 

Institution  of  Mechanical  Engineers,  Vo  I  „  155,  PP .  211~ol  (1946) 

3.  Rusetskiy,  A. A.  and  Favorskaia,  A„V„,  "Sidrodinamicheskie 
Momenty  Nalopasr.lakh  Grebnykh  Vintcv  Reguliruemcgo  Shage,"  Sudostroenie, 
No.  12,  pp.  7-1 2  (1956) 

4.  Van  de  Voorae,  C.B.,  "De  Berejkening  van  het  Bladverstelkoppel 
BiJ  Verstelbare  Schoeven  en  de  Middlen  Ter  Verkleining  van  dit  Koppel," 
Schip  en  Werf,  Vol .  26,  No.  17,  pp.  5C>5- 5^9  (21  August  1959) 

5«  Riegels,  F.,  "Uber  Die  Berechnung  der  Druckverteilung  vcn 
Profilen,'*  Teclmifache  Berichte,  Von.  iG  (1943);  aloe  David  Tay'.rr  Model 
Basin  Report  Aero  Memorandum  28  (March  1955) 

6.  Abbott,  I.H.,  et.  al.,  "Summary  of  Airfoil  Data,"  National 
Advisory  Committee  for  Aeronaut  ics  Report  No,  824  ( 1945  - 


29 


7°  Nonweiler,  T„,  "'The  Design  of  Wing  Sections,"  Aircraft  Engineering, 
pp.  216-227  (July  1956) 

8.  Eckhardt,  M.K.  and  Morgan,  W.B.,  "A  Propeller  Design  Method," 

Society  of  Naval  Architects  and  Marine  Engineers,  Vol.  63,  pp»325“37^> 

(1955) 

9„  "Explanatory  Notes  for  Resistance  and  Propulsion  Data  Sheets," 

The  Society  of  Naval  Architects  and  Marine  Engineers  Technical  and  Research 
Bulletin  No,  1-13  (1952) 

10.  Hecker,  R.,  "Manual  for  Preparing  and  Interpreting  Data  of  Propeller 
Problems  Which  Are  Programmed  for  the  High-Speed  Computers  at  the  David 
Taylor  Model  Basin,"  David  Taylor  Model  Basin  Report  1244  (Aug  1959) 

11.  Morgan,  W.B.,  "An  Approximate  Method  of  Obtaining  Stress  in  a 
Propeller  Blade,"  David  Taylcr  Model  Basin  Report  919  (Oct  1954) 

12.  Burrill,  L.C.,  "Calculation  of  Marine  Propeller  Performance 
Characteristics,"  North  East  Ocas*.  Institute  of  Engineers  and  Shipbuilders, 

pp.  269-293  (1943-1944) 

13.  Kervin,  J.E.,  "Machine  Computation  of  Marine  Propeller  Charac¬ 
teristics,"  International  Shipbuilding  Progress,  Vol.  6,  No.  60  (Aug 
1959) 

14.  Lodkin,  A.S.  and  R^setskiy,  A. A.,  "Towards  the  Calculations  of 
the  Hydrodynamic  Characteristics  of  Varible  Pitch  Propellers,"  Bureau 
of  Ships  Translation,  No.  J26  (196c)  [Translated  from  Sudostroenie, 

Vol.  11,  No.  7,  pp.  9-12,  i960] 


% 


* 


& 


APPENDIX 


DERIVATION  OF  EXPRESSION  FOR  CENTRIFUGAL  COMPONENT  OF  SPINDLE  TORQUE 
The  centrifugal  force  on  an  element  of  blade  mass,  dm,  is  given 
by  (see  Figure  4):  , 

dCF 

dCF  ncol(  JZ)(  dnXad  If) 

The  component  of  this  centrifugal  force  in  the  plane  normal  to  a  radial 
axis  is: 

SIN  X  c\CF  = cjx(  A  2?)(dn.)(n d  y)  a  sim  ^ 
where  X  is  measured  from  the  axis  in  question.  Thus  the  spindle  torque 
due  to  the  centrifugal  force  on  an  element  of  blade  mass  is: 


Q"  =  JCF  SIN  X  z  - £0* n.  SlNX  dn  adY  dZ 
integrating  over  the  sectional  area: 

Q'c  =/°t>  uxAnJfxL  n  5/M  X  Al  ndx 

Evaluation  of  this  integral  over  the  blade  sections  would  be  a  lengthy 
and  tedious  process.  If  s^X  is  replaced  by  X  the  integral  takes  the 
form  of  the  product  of  inertia  of  the  developed  blade  sections  and  thereby 
use  can  be  made  of  the  moments  of  inertia  of  the  sections^  which  greatly 
reduces,  the  amount  of  labor  required  to  calculate  the  centrifugal  com¬ 
ponent  of  blade  spindle  torque.  Setting  X  ~  SI>J  >  Muation  [a] 
becomes: 


R 

Q[  ntr  J2-  ndX 


[b] 


setting  Y  =  a  X  , 

(%-/ohu'-dn.jfY2  dYdi 
But  jj YZ  dY 62  *  Iyl 


where  I  is  the  product  of  inertia  of  the  developed  sectional  area 
about  the  spindle  axis  as  the  origin  of  the  Y  &  Z.  axes.  The  product  of 


inertia  about  the  spindle  axis  is  now  expressed  in  terms  of  the  product 
of  inertia  about  the  section  centroid  as: 


iyi  -  ryJ  +  Firi?  [el 

The  product  of  inertia  about  the  section  centroid  is  next  expressed  in 
terms  of  the  maximum  and  minimum  moments  of  inertia  through  the  section 
centroid: 

T  _  (^MRX  3-Mln)  SIN  Z(j) 

hi  2 -  ■  [fl 


(This  is  under  the  assumption  that  the  major  and  minor  axes  of  the  section 
are  parallel  and  perpendicular,  respectively,  to  the  section  chord.) 
Substituting  Equations  [d],  [e],  and  [f]  into  Equation  [c]  yields  the 
desired  relationship: 


)  sixj 

Z 


RY2 


[10] 


•  The  negative  sign  is  included  so  that  the  sign  convention  will  agree  with 
that  ,of  the  hydrodynamic  component  of  spindle  torque;  that  is,  positive 
spindle  torques  being  those  tending  to  twist  the  blade  towards  larger 
positive  pitch  settings. 

In  order  to  give  an  idea  of  the  magnitude  of  the  error  introduced 
by  the  approximation  ft  =  sixi  %  ,  a  comparison  is  made  of  Equations  [a] 
and  [b]  (exact  and  approximate  solutions,  respectively)  evaluated  over  a 
simple  geometric  section.  Consider  the  rectangular  section: 


Exact  solution  (Equation  [a]): 


U)zlZzdri 


2  sinX  dXdZ 


SIN  X  JX c/Z  -y/k.  (~2~~  Xcos  Xj 


•  Approximate  solution  (Equation  [t>]): 


Qc  =yOb  u)z  nz  cln  JJ  Z  X  JZ  dY 

a  -w  A*) 


Exact  Q'c 
Rpprox  Q'c 


siN‘l)a 


Sxo-ct  Q;  =  smz  §  -  5IM1  $ 
Rpp-xox  Q' 


I  im 

K~>K 


SIN * 


sin* 


sin  y, 

* 


COS  ^  ) 
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Exact  Q' 


Approximate  Q' 


for  rectangular  section  as  function  of  X,  and  ^ 


These  errors  are,  of  course,  larger  than  would  occur  for  a  blade  section 
with  similar  ^  and  >  since  the  thickness  of  a  usual  blade  section 

is  reduced  where  X  is  largest . 

On  the  propeller  used  in  the  example  calculation  of  this  report,  a 
maximum  X  of  40°  occurs  (at  the  section  at  the  hub).  At  greater  radii 
where  the  centrifugal  force  is  larger,  the  maximum  sectional  X's  are 
smaller.  This  indicates  that  the  approximation  %  =  siw  X  introduces  an 
error  of  less  than  5$  in  the  centrifugal  component  of  spindle  torque 
for  this  propeller.  The  value  of  X  on  controllable-pitch  propellers  as 
a  rule  is  not  exceedingly  large,  since  the  blade  width  is  limited  by 
interferences  between  the  blades  at  various  pitch  settings  and  greatly 
skewed  blades  are  not  generally  used.  For  any  particular  propeller,  the 
reader  may  decide  whether  or  not  the  approximation  is  close  enough  for  hi's 
purpose.  If  not,  graphical  or  numerical  integration  of  Equation  [a]  may 
be  used  to  obtain  a  more  accurate  solution. 
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